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Abstract 

A general analysis of Q-ball solutions of the super symmetric F-term hybrid 
inflation field equations is given. The solutions consist of a complex inflaton field 
and a real symmetry breaking field, with a conserved global charge associated 
with the inflaton. It is shown that the Q-ball solutions for any value of the 
superpotential coupling, k, may be obtained from those with n = 1 by rescaling 
the space coordinates. The complete range of Q-ball solutions for the case k = 1 
is given, from which all possible F-term inflation Q-balls can be obtained. The 
possible role of F-term inflation Q-balls in cosmology is discussed. 
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1 Introduction 



Hybrid inflation models are widely regarded as promising candidates for inflation in 
the context of supersymmetry (SUSY). There are two classes of SUSY hybrid inflation 
model \X\, F-term [2] and D-term inflation depending on which part of the scalar 
potential drives inflation. An important feature of all inflation scenarios is the process 
by which inflation ends and the Universe reheats. In SUSY hybrid inflation models 
this is likely to occur via rapid tachyonic growth of sub-horizon quantum scalar field 
fluctuations into a space-dependent classical field, a process known as tachyonic pre- 
heating jlj. The subsequent evolution of the resulting inhomogeneous classical field 
is a complex numerical problem. In a recent simulation of the evolution of the scalar 
fields at the end of SUSY hybrid inflation [S] , the formation of oscillons [H] Ej made 
of the inflaton sector fields (essentially droplets of Bose condensate held together by 
an attractive interaction) was observed. The possibility of the formation of such non- 
topological solitons at the end of SUSY hybrid inflation was first suggested in jH] (where 
they were called 'inflation condensate lumps') and their formation was also observed in 
a numerical simulation of the growth of perturbations of the symmetry-breaking field 
in SUSY hybrid inflation models Similar oscillon states have been observed in 
numerical simulations of the growth of scalar field fluctuations along SUSY flat direc- 
tions JU] and in a chaotic inflation model JT] based on the SUSY flat direction scalar 
potential [12] . In those cases the oscillons were observed to subsequently decay to Q- 
ball anti-Q-ball pairs. This is a natural possibility since the Q-ball is an energetically 
preferred stable state, to which the unstable neutral oscillons will naturally evolve via 
the growth of perturbations of the phase of the complex scalar field. A similar decay 
of hybrid inflaton oscillons to Q-balls might therefore be expected if stable Q-balls 
made of inflaton sector fields exist, resulting in a Q-ball dominated post- inflation era. 
This would have important consequences for reheating [HJ EI] and post- inflation scalar 
field dynamics |14j . 

In a previous paper we demonstrated the existence of Q-ball solutions of the field 
equations of the D-term inflation model [T^j. It was noted that, for a particular choice 
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of couplings, the Q-ball solutions in the D-term inflation model are equivalent to Q- 
balls in the minimal F-term inflation model, based on a single symmetry breaking 
field which carries no gauge charge. On the other hand, the most commonly studied 
F-term inflation model is that based on a pair of symmetry-breaking fields, $ and 
$, transforming as a conjugate pair under a gauge group. (In the following we will 
refer to this as standard F-term inflation.) The Q-ball solutions in this model are not 
equivalent to those in D-term inflation. 

In this paper we present a detailed analysis of Q-ball solutions in the standard 
F-term inflation model. We will show that the Q-ball solutions for any superpotential 
coupling, re, can be made equivalent to those with re = 1 by rescaling the space 
coordinates. As a result, we are able to study the properties of F-term inflation Q-balls 
in a completely general way. We will focus on the case of global SUSY throughout. The 
effect of supergravity (SUGRA) corrections in realistic models will be briefly discussed 
in our conclusions. 

The paper is organised as follows. In Section 2 we derive the equations which 
minimize the energy for a fixed global charge ('Q-ball equations') for the standard 
F-term inflation model. In Section 3 we show that the Q-ball solutions for any super- 
potential coupling can be obtained by rescaling the space coordinates and solving the 
field equations for the case re = 1. In Section 4 we present numerical solutions of the 
standard F-term inflation Q-ball equations and give the full range of solutions for the 
case re = 1. In Section 5 we present our conclusions and discuss the possible role of 
F-term inflation Q-balls in cosmology. 



2 Q-Ball Equations of the Standard F-term Infla- 
tion Model 



The most commonly studied form of F-term inflation model has the superpotential 
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where S is the inflaton and $ and $ are the symmetry-breaking fields. We may choose 
fi 2 to be real and positive. The symmetry-breaking fields typically represent oppositely 
charged gauge multiplets in realistic models. The scalar potential is then 
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where $ and $ now represent the components of the gauge multiplets which gain an 
expectation value. In this we have included a generic D-term for the case where $, 
$ transform under a gauge group. The vacuum state will correspond to S = and 



$1. 



= fi 2 . In the case where there is a D-term, the vacuum will be such that 
We can choose $ to be real and positive, in which case the vacuum expectation values 
(VEVs) will correspond to $ = $ = [l. We will focus on this case in the following. 

A superpotential proportional to S, together with a minimal Kahler potential term, 
is essential in order to avoid the ^-problem in F-term inflation models once SUGRA 
corrections are included [16j. A superpotential which is proportional to S implies that 
there is a global U(l) R-symmetry, under which only the inflaton transforms, and an 
associated global charge. The existence of a conserved inflaton charge, Qs, implies the 
possibility of Q-ball solutions of the F-term inflation field equations. 

For a model with a single complex scalar field, and a global U(l) symmetry under 
which \1/ — > e ta ty, the Q-ball configuration is derived by minimizing the functional 



E u {V,V,u) = E + u(Q- / d 3 xp Q 
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where uo is a Lagrange multiplier fixing the charge, E is the total energy of the field 
configuration 
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and pQ is the global charge density 

p Q = iffity - \I/tvj/) . 
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This minimizes the energy for a fixed global charge, Q. E w may be equivalently written 

as 



\^-ilu^\ 2 + |V*r + V(V) -lu 2 \^\ 2 d 6 x + uQ . (6) 



Minimizing with respect to ^ implies that \P(x, t) = \l/(x)e lw *. Substituting this into 
Eq. (0) gives 



E u (*(x),u) = J [|W(x)| 2 + ^(x)) -cu 2 |^(x)| 2 ] d 3 x + toQ . (7) 



Extremizing this with respect to ^(x) then implies that 
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where = V(^) — a; 2 |\I/| 2 . ^(x) could still have a space-dependent complex phase, 
#(x). However, if V(^f) = V(\^f\) (as it must when \1/ transforms under a £7(1) sym- 
metry), then Eq. (|SJ) is generally minimized by the choice 9 = constant, which may be 
chosen such that ^(x) is real. A minimum energy configuration should be spherically 



We will refer to this as the Q-ball equation. The solutions of Eq. © should satisfy the 
boundary conditions that the field tends to the vacuum as r — > oo and that dip/dr — > 
as r — > 0. On physical grounds we expect that there will be a unique Q-ball solution 
of this form for a given global charge. The resulting solution will be a stable Q-ball 
if E/Q < my, since it would then require additional energy to remove the conserved 
global charge into the vacuum. 

The procedure for solving Eq. © numerically to obtain Q-ball solutions for a given 
value of uj is as follows. We first input a value for ip(r — 0) = ip Q , impose the boundary 
condition that dip/dr = at r = 0, and solve the Q-ball equation for ip(r). We then 
vary ip until we obtain a solution which asymptotically approaches the vacuum as 
r — > oo. This will correspond to a possible Q-ball solution. Finally, we calculate E 
and Q for this solution and check that E/Q < my. In practice we can only check 
the asymptotic approach to the vacuum out to a finite radius, r en d. So long as r en d is 
large compared with the dynamical scale of the scalar field, m^ 1 , this will give a good 
approximation to the true Q-ball solution. Larger values of r end require progressively 
more accurate values of i/) . 



symmetric. Therefore, with \l/(x) = ip(r)/\/2 (where ip(r) is real), Eq. (jHJ) becomes 

d 2 ijj 2 dijj dV 2 
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The above analysis generalises to the case of SUSY hybrid inflation with two or 
more complex scalar fields. The Q-ball equations generally become a system of coupled 
equations for the inflaton, s, and symmetry-breaking field, 0, with a single Lagrange 
multiplier associated with the inflaton. The numerical procedure is as before, except 
that we must scan through values of s D and O , looking for solutions which satisfy the 
requirement that s and <ft asymptotically approach their VEV as r — > oo. 

For the case of standard F-term inflation there are three complex scalar fields. 
However, we can only associate a Lagrange multiplier with fields which have vacuum 
expectation value equal to zero. This follows since the energy minimization with a 
Lagrange mulitplier 7 for a field $ which transforms under a global U(l) symmetry 
implies that the Q-ball solution must be of the form $(x, t) = 0(r)e n */v / 2- This is 
only compatible with a constant VEV as r — > 00 if <p(r) — > 0. Therefore in the case 
of F-term inflation, where $ and $ have non-zero VEVs, we can only associate a 
Lagrange multiplier with the S field. 

Thus we will look for Q-ball solutions with a global charge Qs associated with the 
R-symmetry under which S — > e ia S. The energy functional is then 
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where 

^(|S|,$,$) = n|S|,$,$)-^ 2 |S| 2 , (11) 

with V(|5|, $) given by Eq. (J2J). Minimizing E w with respect to the time derivatives 
then implies that S(x,t) = S(x.)e iujt , $(x,£) = $(x) and $(x,t) = $(x). Eq. (JE) then 
becomes 

E u = f [|V5(x)| 2 + |V$(x)| 2 + |V$(x)| 2 + K,(|5|,$,$)] d 3 x + uQ s . (12) 

The potential is explicitly dependent upon the phase of $ and $. However, by inspec- 
tion of Eq. 0, the minimum of E w with /i 2 real and positive will generally correspond 
to $ and $ real and positive. Therefore with S'(x) = s(r)/\/2, ^*( x ) = 0( r )/"\/2, 
$(x) = 0(r)/-y/2, the Q-ball equations analogous to Eq. Q from extremizing Eq. (JT2j) 
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Since in the case with a D-term the VEV of $ and $ are both equal to //, we will look 
for Q-ball solutions with 0(r) = 0(r). In this case the equations reduce to 

d 2 s 2 ds 



2 j 2 2 
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The energy and charge of the resulting Q-ball soultion are given by 
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where 
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and 



= J PQ s d 3 x = oj J Airr 2 s 2 dr , 



(20) 



where pq s = i (S^S — S^S^j . There is no factor of 1/2 in front of (d(f)/dr) 2 contribution 
to E, since this includes the energy associated with in the case where = 0. 



3 Scaling Property of the Q-Ball Equations 

The two free parameters of the standard F-term inflation superpotential are the mass 
scale, /x, and the superpotential coupling, ft. The mass scale may be eliminated by 
using units such that /i — 1. We can avoid having to produce a set of Q-ball solutions 
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for each value of re by rewriting the Q-ball equations in terms of r = nr. The Q-ball 
equatons then become 



d 2 s 2 ds 
df 2 f df 
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where Q — ou/k. These equations have the same form as the Q-ball equations for the 
case re — 1 in /i — 1 units. The energy and charge are given by 
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where 



and 



Qs = %l ^r 2 s 2 df ee % . 



K 



K 



(24) 



(25) 



E and Qs are the energy and charge of the corresponding k — 1 Q-ball. The inflaton 
mass is given by m s = V2k,/jl. Therefore the Q-ball stability condition, E/Q s m s < 1, 
is independent of re. 

Since all the Q-balls for an arbitrary value of re have a one-to-one mapping to those 
with re = 1, the full set of Q-ball solutions for the case re = 1 will constitute a complete 
solution for the range of possible Q-balls in standard F-term inflation. 



4 Numerical Q-ball Solutions 

In Table 1 we give the set of re = 1 Q-ball solutions we have generated in order 
to analyze the general properties of F-term inflation Q-balls. We use fi — 1 units 
throughout. The inflaton mass is then m s = V2k and the VEV of equals 
The solutions were generated using a second-order Runge-Kutta routine with radial 
increment Ar = 0.01. 
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UJ 


So 


4>o 


E 


Qs 


E/Q s m s 


R 


^ G-Tld 


0.58 


7.7891 


1.03 x 10~ 6 


2756.42 


3654.07 


0.533 


5.06 


6.40 


0.6 


7.2700 


5.40 x lO" 6 


2370.57 


3006.14 


0.557 


4.92 


8.0 


0.65 


6.6699 


5.24 x 10~ 5 


1855.8 


2180.4 


0.602 


4.57 


8.0 


0.7 


6.1518 


3.16 x 10~ 4 


1479.41 


1621.36 


0.645 


4.28 


8.0 


0.8 


5.2854 


4.1436 x 10~ 3 


980.86 


951.66 


0.729 


3.83 


8.0 


0.9 


4.57364 


2.3758 x 10~ 2 


681.07 


596.90 


0.807 


3.51 


8.0 


1.0 


3.95008 


0.08160 


491.15 


396.03 


0.877 


3.32 


8.0 


1.1 


3.35932 


0.202129 


366.43 


276.78 


0.936 


3.24 


8.0 


1.17 


2.93500 


0.334867 


304.56 


222.17 


0.969 


3.28 


8.0 


1.25 


2.38244 


0.552551 


255.58 


181.67 


0.995 


3.61 


6.90 


1.30 


2.0631 


0.69888 


222.45 


156.02 


1.013 


3.60 


6.30 



Table 1: Properties of F-term inflation Q-ball solutions with k = 1. 



In Figure 1 we show s(r) and </>(r) for a typical k = 1 F-term inflation Q-ball. 
In this case (uj = 1.1) we find that s(r) can be approximately fitted by a Gaussian, 
s(r) = s G exp(— r 2 /R 2 ) with s Q = 3.36 and R = 2.24. This is a good fit for r < 4.0, but 
for larger r the Gaussian gives a value of s(r) somewhat smaller than the numerical 
solution (by a factor of 10 at r = 6.0). Nevertheless, the Gaussian provides a useful 
analytical expression for s(r), which can be used when studying the cosmology of 
Q-balls HI]. 

In Figure 2 we show the ratio E/Qsms as a function of uj. As discussed below, there 
is a lower limit on u, approximately 0.57, at which value the energy and charge of the 
Q-ball diverges. An upper limit on uj is imposed by the requirement that E/Qsms < 1 
for a stable Q-ball. From Figure 2 we see that this requires that to ^ 1.26. (We will 
apply this upper limit to the remaining figures.) 

In Figure 3 we show the energy and charge of the Q-ball solutions as a function 
of uj. Both E and Qs display vertical asymptotic behaviour as uj approaches 0.57, 
therefore Q-balls of arbitrarily large energy and charge are possible. The upper bound 



8 



0123456789 

r 

Figure 1: F-term Inflation Q-ball for the case k = 1, u = 1.1. 

on uj coming from Q-ball stability imposes a lower bound on the energy and charge 
for a stable Q-ball solution to exist, E min pa 250 and Q m in ~ 180. 

In Figure 4 we show O as a function of uj. We find that, as uj decreases, the value 
of <p becomes small and asymptotically approaches zero as uj —>■ 0.57. The extremely 
small and precise values of <p required as u approaches 0.57 makes it difficult to find Q- 
ball solutions in this limit. The value of s Q also begins to increase rapidly as uj — > 0.57, 
as shown in Figure 5. 

In Figure 6 we show the Q-ball radius, R, as a function of uj, where we have defined 
R to be the radius within which 90% of its energy resides. We find that the radius 
of the Q-balls is relatively insensitive to the energy and charge, varying from 3.27 to 
5.06 over the entire range of Q-ball solutions. 

For the case of Q-balls with k < 1, the effect of rescaling the k = 1 coordinates 
will be to increase the lower bound on the energy to E min pa 250/k and the charge 
to Qmin pa 180/k 2 . The corresponding Q-balls will have a larger physical radius, 
R = R(k = The range of u for which stable Q-ball solutions exist becomes 

0.57k < u < 1.26k. 
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Figure 4: <\> Q versus uo for k = 1 Q-balls 
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Figure 5: s Q versus w for k = 1 Q-balls 
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Figure 6: Radius versus u for k = 1 Q-balls 

5 Conclusions and Discussion 

In this paper we have studied the complete range of possible Q-ball solutions in the 
standard F-term inflation model. The solutions are parameterized by the frequency of 
rotation of the inflaton in the complex plane, u. It is useful to restate our main results 
in terms of the physical inflaton mass (ms = \/2k in /x = 1 units). For a superpotential 
coupling k, stable Q-balls exist for 0.40ms < u ^ 0.89m,g. As uo — > OAOms, the charge 
and energy of the Q-ball diverge. There is a lower bound on the charge and energy 
for which a stable Q-ball exists, Q m in ~ 180/k 2 and E min > 180ms/ k 2 . The radius of 
the Q-ball, within which 90% of its energy resides, is between 4.62m^ 1 and 7.16ms 1 . 

The existence of Q-ball solutions of the standard F-term inflation model, combined 
with the rapid growth of quantum scalar field fluctuations at the end of inflation, 
implies that Q-ball formation at the end of F-term inflation is a possibility. What 
can we deduce about F-term inflation Q-ball cosmology, should their formation at the 
end of inflation be confirmed in future numerical simulations? The typical energy and 
charge of such Q-balls, assuming that uo is not close to its upper or lower limit, would 
be E m 10 3 ms/K 2 and Qs ~ 10 3 /k 2 . It was recently shown that cosmic microwave 
background constraints impose limits on the range of K in F-term inflation models, 
10~ 7 ^ K ^ 10~ 2 [17J. Within this range, the global charge carried by the Q-balls can 
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be very large, up to Qs ~ 10 . In the case where gravity-mediated SUSY breaking 
A-terms exist during inflation, the lower bound becomes k ^ 10 -5 |18j . corresponding 
to Q-ball charges up to Qs ~ 10 13 . Such large charge Q-balls would be expected to 
have interesting consequences for reheating, delaying the process relative to the case 
of a conventional perturbatively decaying homogeneous inflaton condensate [T21 119j. 
The consequences of inhomogeneous reheating occuring via Q-ball decay would be 
particularly interesting if it occurs close to a physically significant temperature, such 
as that of the electroweak phase transition or dark matter freeze-out. 

The existence of a lower bound on the charge of a stable Q-ball has interesting 
implications for the Q-ball decay process. The inflatons making up the Q-ball would 
decay perturbatively until the charge drops below Q m i n , at which point the Q-ball 
would rapidly expand due to its gradient energy and dissociate into inflatons. 

It should be emphasized that the existence of Q-ball solutions is not dependent 
upon the existence of a symmetry-breaking phase transition in hybrid inflation models. 
For example, it is straightforward to show that Q-ball solutions also exist in smooth 
hybrid inflation models [20], where inflation ends without any phase transition |21j . 

In this paper we have considered F-term inflation Q-balls in the limit of global 
SUSY. The F-term inflation Q-ball solutions are stable as a result of a global U(l) 
symmetry, which must correspond to an R-symmetry in order to be consistent with the 
K[i 2 S term in the superpotential. However, in realistic models we must also consider 
the effect of SUGRA corrections. If SUGRA with hidden sector R-symmetry break- 
ing is considered then there will be U(l)-breaking SUGRA corrections which could 
destabilise the Q-balls 1 . For example, in the case of gravity-mediated SUSY breaking, 
an R-symmetry breaking A-term of the form Va ~ m^^K^S + h.c. will generally be 
added to the scalar potential jTH] , which would add a constant term to the right hand 
side of the S field equation (Eq. (fTBJO of the form m 3 / 2 K/i 2 . The magnitude of s and 
4> inside the Q-ball is typically of the order of n, so the ratio of the U(l)-breaking 
A-term to the U(l)-symmetric terms in Eq. (|16|) is approximately m^n/ms- Since 
1 We note that in the case of D-term inflation the Q-balls JS] are stable due to a global U(l) 
symmetry of the superpotential and so are not affected by R-symmetry breaking SUGRA corrrections. 
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the inflaton mass, mg, is typically many orders of magnitude larger than m.3/2, the 
A-term will make a very small contribution to the S field equation compared with the 
U(l)-symmetric terms. In this case there will be a time- dependent oscillon solution 
with properties very close to those of a Q-ball. Given the small size of the A-term 
and the inherent stability of oscillon solutions jHj, we would expect the lifetime of 
the Q-ball-like objects to be very long compared with the dynamical timescale frig 1 . 
However, a numerical solution for the effect of the A-term on the Q-ball lifetime will 
be necessary in order to find out whether it could dominate over perturbative inflaton 
decay. 

Beyond cosmology, the F-term inflation Q-ball provides an interesting example of 
a non-topological soliton in the context of a physically-motivated scalar field theory. 
It would be interesting to try to understand analytically some of the features which 
we have found numerically, such as the lower bound on the charge and energy and the 
finite range of u. The properties of the Q-balls might also shed light on the dynamics 
of the related oscillons which occur in the case of a real inflaton field. For example, the 
existence of a lower bound on the Q-ball energy to form a stable Q-ball might imply 
a similar lower bound on the oscillon energy. The advantage of the Q-ball is that it 
is possible solve exactly for its time-dependence, whereas oscillons generally have a 
complicated time- dependence which is difficult to study numerically 
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